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Some of the assumptions used in treatment of the mode theory of VLF radio propagation 
are discussed briefly. The connections with geometrical-optical theories are also pointed out. 



Iii several recent papers, the author has discussed 
waveguide theory of radio propagation for a smooth 
spherical earth and a concentric anisotropic iono- 
sphere [Wait, I960; L963a, L963b]. To facilitate the 
solution of this intrinsically difficult problem, certain 
assumptions were made. In view of great interest 
in this problem, it appears to be worthwhile to 
discuss briefly the nature of these assumptions since 
they have been alluded to in a series of three papers 
[Berry, 1963; Jollier and Berry, 1963; Johler, 1963]. 

The source is taken to be a vertical electric dipole 
located at 6=0 and r = b of a spherical coordinate 
system. The earth's surface is located at r=a and 
the reference reflecting layer is at r=c. The 
medium between a and c is regarded as free space. 

Within certain limitations, which have been 
discussed extensively [Wait, 1962], the boundary 
conditions at r=a and r=c may be described in 
terms of surface impedances. In the general case, 
the surface impedances are functions of the (complex) 
angle of incidence. However, for treating VLF 
propagation in the earth-ionosphere waveguide, it is 
justified to regard the elements of the surface im- 
pedance matrices as constants. As a result, lateral- 
type waves are neglected. The latter correspond to 
waveguide modes propagating within the earth and 
the ionosphere. In the limiting case of a flat earth 
model, they become part of the continuous spectrum, 
and they are usually described as the branch-cut 
waves. A detailed investigation of these lateral 
waves for a planar model has been carried out by 
Anderson [1962]. As expected [Wait, 1960], they 
do not contribute significantly to the total field. 

Another consequence of using impedance boundary 
conditions is that waves "reflected" at the center of 
the earth are neglected. In essence, this is ac- 
complished by Bremmer [1949, 1958] when he re- 
places the spherical Bessel functions j n (k g r) by 
hn } (Jc g r) , where k g is the propagation constant of 
the earth. In view of the large magnitude of the 
imaginary part of k g r, it turns out that this is an 
excellent approximation. 1 



1 For purposes of this discussion, the time factor is exp (forf). 



Thus, in one fell swoop, the hit (Mai waves and the 
waves reflected from the center of the earth are 
dispensed with by the introduction of impedance 
boundary conditions. Furthermore, in the applica- 
tion of the Watson transform, the desired even 
property of the integrand is assured [Wait, 1960, 
1962]. 

In much of the theoretical development of this 
subject, the spherical wave functions h^ l) (x) and 
hi, 2) lx), of genera] argument x, are ultimately approxi- 
mated in terms of Airy functions (or Hankel func- 
tions of one third). Thus, if \v — x|<<V /3 , it is 
permissible to write 



and 
where 



xh { v l) (x)^-i(x/2) l / & w 2 (t), 

xhi 2) (x)^+i(x/2) 1/6 w l (t), 
t=(v+h-x)(2/xy<\ 



and w 2 (t) and w Y (t) are Airy functions. For example 



W : 



2 (t)=^[Bi(t)+iAi(t)], 



wi(t) =\Tr[Bi(t) —iAi(t)], 

in terms of Airy functions Ai(t) and Bi(t) defined 
by Miller [1946]. 

These Airy function representations are highly 
accurate for the VLF mode problem as has been 
demonstrated by Wait [1960] and Spies and Wait 
[1961]. Using such approximations, the modal 
equation takes on a reasonably tractable form, and 
numerical results for attenuation and phase velocity 
may be computed in an economical manner. Despite 
this, some further approximations are often war- 
ranted if simplicity and physical insight are to be 
gained. For this reason, the Debye approximations 
for the spherical wave functions are found to be 
very useful [Wait, 1962]. In the present context, 
they correspond to utilizing the leading term in the 
asymptotic expansion of the Airy functions valid for 
large negative t. 

The Debye approximated mode equation is 
accurate when \C n (ka/2)*\ and |(7;(&a/2)*|>2, where 
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C n and C' n are the cosine of the (complex) angle of 
incidence for the earth and the ionosphere, respec- 
tively. Because (ka/2)^, at VLF, is the order of 
20 and C' n is never smaller than about 0.15, the 
latter condition is never seriously violated. How- 
ever, the quantity \C n {kaj2)*\ may become very 
small at VLF for important modes at frequencies 
greater than 10 kc/s. 

Another approximation occasionally used corre- 
sponds to the assumption of a flat earth. Such a 
model has been used extensively by Budden [1962]. 
The resulting simplicity permits a careful study of 
the nature and the properties of the modes. Actu- 
ally, the flat-earth mode equation is obtained as 
a limiting case of the spherical-earth equations when 



both |tf„(£a/2)*| and 



/12 



>>1. Extensive com- 



parisons between flat-earth and spherical-earth 
computations have been reported previously [Wait 
and Spies, 1960; Wait, 1962]. 

In certain cases, when the convergence of the sum 
of waveguide modes is poor, an alternative repre- 
sentation is desirable. Such a representation may 
be found from geometrical optics wherein the total 
field is considered to be a sum of rays or hops. The 
ray or hop series has been obtained directly from 
the contour-integral representation by Bremmer 
[1949, 1958]. This result is derived by expanding 
the integrand, which is of the form [1— R g R i ]~ 1 } into 
a geometrical progression. Individual terms are 
then evaluated asymptotically to yield the classical 
geometrical-optical formulas. The method has been 
extended by Wait [1961] to permit the results 
to be used in the vicinity of the geometrical horizon 
of the various rays (i.e., hops). A key point in the 
development is the ordering of terms so that R g and 
Ri are defined as (spherical) reflection coefficients 
connecting "downgoing" and "upgoing" wave func- 
tions. The latter have the forms h^(kr) and 
hi 2) (kr), respectively, in the earth-ionosphere region. 
If another choice of wave functions is used (e.g., 
j v (kr) and hj, 1} (kr)), the individual integrals in the 
expansion should not be identified with geometrical- 
optical hops in the usual sense. The distinction 
between these methods of expansion is not trivial. 
In the first method, the only significant poles are 
those of the groundwave type, whereas, in the 
second method, additional poles appear which are 
closely related to the usual waveguide modes. Thus, 



if care is not exercised in the choice of wave functions, 
the individual integrals in the expansion will corre- 
spond to hybrid waves which have both guided-wave 
and ground-wave components. While the resultant 
field as computed should be identical, the utility of 
such hybrid expansions could be questioned. 

For further discussion of these various points, the 
reader might wish to consult the referenced literature. 
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